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We argue a scenario motivated by the context of string landscape, where our universe is produced 
by a new vacuum bubble embedded in an old bubble and these bubble universes have not only 
different cosmological constants, but also their own different gravitational constants. We study 
these effects on the primordial curvature perturbations. In order to construct a model of varying 
gravitational constants, we use the Jordan-Brans-Dicke (JBD) theory where different expectation 
values of scalar fields produce difference of constants. In this system, we investigate the nucleation 
of bubble universe and dynamics of the wall separating two spacetimes. In particular, the primordial 
curvature perturbation on superhorizon scales can be affected by the wall trajectory as the boundary 
effect. We show the effect of gravitational constant in the exterior bubble universe can provide a 
peak like a bump feature at a large scale in a modulation of power spectrum. 

PACS numbers: 98.80.-k, 98.90.Cq 


I. INTRODUCTION 

In the context of the string theory landscape [ij, 
there are various possibilities that multiple vacua can 
be produced whose semiclassical tunnelling process from 
metastable vacua to new vacua can occur. In this re¬ 
spect, it can allow us to argue that nucleation of spher¬ 
ical symmetric regions, namely Bubble in space which 
are constructed by a new vacuum and expand into the 
old vacuum. One of any bubble universe created by such 
first-order phase transition may be regarded as our uni¬ 
verse |2 hI 4|. In particular, inspired by the bubble uni¬ 
verse scenario, an accelerated expansion of the universe 
generated by its vacuum energy would give a natural ex¬ 
planation for inflation realization in the early universe 
as discussed old inflation Q or open inflation SiEl 
scenario in the past. If the bubble universe would be 
inside the old bubble universe, the so-called parent uni¬ 
verse, we are mainly interested in how much information 
we know or detect through any observational signature 
of the parent universe. 

In the viewpoint of Anthropic principle [3, Ell, it is 
possible that several universes take their own choices of 
natural constants. It can allow us to argue a possibil¬ 
ity that landscape vacua have different values of gravi¬ 
tational constants as well as cosmological constants. In 
this paper, the bubble universes have their own different 
gravitational constants and it can give us explore how a 
different gravitational constant in the exterior universe 
affect and relates to physical (especially observational) 
quantities for the observer in the interior bubble universe 
(see also works El for the aim). 

In order to such model building that each bubble uni- 
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verse has its own different gravitational constant, it is 
useful to employ the JBD theory El (see book El for 
review). This theory can achieve a varying gravitational 
constant by identified a vacuum expectation value of the 
scalar field with gravitational constants in the original, 
the so-called Jordan frame. Therefore we will use this 
theory to construct our model and investigate the conse¬ 
quence of the model to our (child) bubble universe, where 
the child bubble is inside the parent bubble with different 
constants. 

Inflationary scenario can give a remarkably successful 
way of generating a consistent spectrum of curvature per¬ 
turbations with recent cosmological observations, such 
as the cosmic microwave background (CMB) a nd g alaxy 
survey to explore a large-scale structure SHI El]. The 
quantum fluctuations generated in a exponential expan¬ 
sion of the early stage of the universe were spread out 
a cosmological horizon scale and their wavelengths are 
then pushed to superhorizon size, which can be reduced 
to classical perturbations seeding cosmic structure. Af¬ 
ter inflation ends, they reenter the horizon in a radiation 
dominated phase and can be detected by such as CMB 
observation (see for review). The recent data by 
Planck satellite gave a perfectly agreement with the 

prediction of inflationary theory. 

In order to discuss an observational consequence of 
bubble universes with different gravitational constants, 
we analyze density perturbations during an inflationary 
phase based on the JBD theory. Note that although 
several inflationary models with/without bubble forma¬ 
tions and the accompanying density perturbations are 
discussed based on the scalar-tensor gravity theories in¬ 
cluding the JBD theory EB“E||, the JBD scalar field is 
assumed to be an inflaton in those models. In the present 
paper, however, we assume that there exists an inflaton 
field in addition to the JBD scalar field, which is used to 
fix only two different gravitational constants. 

If the superhorizon perturbations propagate even 
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though their amplitudes freeze out and approach to the 
boundary of the bubble universe, they would be scattered 
by the layer called bubble wall, separating a new and old 
vacuum. Namely when one join the interior and exterior 
perturbation mode over the bubble wall, a reflected wave 
can be produced. Hence if there would be a possibility 
that one can get information of outside bubble, that is 
only way to analyze perturbations coming back from su¬ 
perhorizon scales outside bubble as a trace of the outer 
region (see [ 2 ^ in a different setup, but in which a similar 
effect has been studied). 

The aim of this paper is to obtain physical signature 
on a primordial curvature perturbation at superhorizon 
scales, which is affected by the boundary of bubble. In 
particular, we can provide how superhorizon mode is 
modified by the effect of different gravitational constant 
in the exterior universe by employing the JBD theory. 

The paper is organized as follows. In section 2 we 
show our model setup and discuss a nucleation of bubble 
universe and a static bubble wall solution. In section 
3, we study dynamics of bubble wall after a nucleation. 
In section 4, we adopt linear perturbation theory to the 
system and match perturbation modes in both sides over 
the bubble wall. In particular, our calculation is the case 
for de Sitter background. Then we discuss the effect of 
different gravitational constant on the obtained results in 
section 5. Section 6 is devoted to conclusion. 


II. MODEL 


We consider the model of bubble universes which have 
their own gravitational constants. The wall of the bubble 
can divide two bubble universes whose constants have dif¬ 
ferent values. We investigate the effect of different gravi¬ 
tational constant on the superhorizon scale perturbation. 

We will introduce the JBD theory [l^, [13] with a po¬ 
tential V(4>) in the Jordan frame as the model of varying 
of gravitational constant, which is expressed as the action 
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where $ is the scalar field whose potential is considered to 
provide their semiclassical tunnelling, namely it has two 
minima called a true and false vacuum respectively. We 
use the units of = 1, where Gjv is the Newtonian 

gravitational constant. 

From such a potential, the bubble universe will be cre¬ 
ated as a standard instanton process. The gravitational 
constant G is related to the vacuum expectation value 
(VEV) of the scalar field as follows 
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The scalar field has a different value depending on each 
potential minima, and then it can read different values 
of gravitational constants. 


In order to analyze this system, it is useful to study dy¬ 
namics in the Einstein frame by using conformal trans¬ 
formation as 5 ^ 1 / = After transformation, 

the action in the Einstein frame is given by 





(2.3) 


Here we introduced a new scalar field: 4> and its potential, 
which are related to the original ones in the Jordan frame 
as 


^ = ^o + y71n$, H(</.) = ^, (2.4) 

where 4>o corresponds to the VEV of the JBD scalar field 
in our universe (<I>o = I). Hereafter we will omit the tilde 
which represents the variables in the Einstein frame for 
brevity. In this paper, as expected to generate a bubble, 
we consider the potential which should give tunnelling, 
which is expressed by 

^(<^) = ^ ■ (2.5) 

Under the approximation that e is small, two vacua are 
obtained by (j)± = ±fj, — where the suffixes 

± describe two minima. In the case where e > 0, the 
tunnelling occurs from 0+ to vacuum. Each potential 
energy is given by V{(j)+) = A -|- G(e^) and U((/)_) = 
A — e -I- G(e^), respectively. 

Since we assume that our universe was created by tun¬ 
nelling from de Sitter universe, our vacuum state is given 
by <(>o = 4 ’-- From the above potential (12.51) . the poten¬ 
tial V($) in the original Jordan frame is found by (12.41) 
. In Fig.l, we plot two potentials V(<I)) and From 

the figure, it is clear that the tunnelling process can be 
realized from a false ( -I- ) vacuum to a true(—) vacuum 


A. Nucleation of Bubble Universe 

We consider a spatially flat Freidmann-Lemaitre- 
Robertson-Walker (FLRW) spacetime in both sides of 
the wall. In this case, the instanton process is calculated 
in the Euclidean spacetime with 0(3) symmetry. In the 
Lorentzian configuration, the line element takes the form 

ds^ = a^{ri){—dri^ + dr^ + r^dil^), ( 2 . 6 ) 

where 77 is a conformal time which is related to a cosmic 
time t hy dt = ad?]. 


^ Note that the possibility for tunnelling is found by use of the 
potential in the Einstein frame. We may misunderstand the pos¬ 
sibility if we look only at the potential in the Jordan frame. See 
an example in the paper 1261 
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FIG. 1: (Left) We plot the potential V($) in the Jordan frame. The tunnelling occurs from $+ to (Center) We plot 
the potential V{(f)) in the Einstein frame. The tunnelling occurs from <()+ to (f)-. (Right) We plot the symmetric potential 
(12.111) and asymmetric one (12.131) . The positive vacuum 0+ where the local minimum is located moves to the right direction 
by 6-times with the same height of the potential barrier. 


Following 0,0, [il, the calculation has been done by 
use of the complex time path approach, in which dynam¬ 
ics is studied by the complex time rj. By adopting the 
thin-wall approximation, in which a thickness of wall is 
small enough compared to the size of a bubble, the bub¬ 
ble radius r(? 7 ) is solved as a function of rj. The action is 
evaluated as 


S = 


dr] 


-ATraa^{r])y{r])^Jl- {dr,r{r]))^) . (2.7) 


This action contains both contributions from the volume 
and surface area of three-dimensional sphere. Here we 
have used e and cr, which are the energy density differ¬ 
ence between two vacua and the surface energy density 
of the wall, respectively. The equation of motion is given 
by taking a variation of the action S with respect to r. 
Solving it, we find a classical trajectory of the wall r{T]). 
In the backward of time rj, we trace a shrinking bubble 
radius and reach to a turning point where the canoni¬ 
cal momentum p = dC/d{drir) vanishes. At this turning 
point (ry = iji), a bubble nucleation occurs. By the an¬ 
alytic continuation, we find the wall trajectory in the 
complex T] plane, which shrinks smoothly to zero size of 
a bubble. In order to perform this procedure, it is useful 
to rewrite the equation as one for r]{r). The task to do 
is reduced to solve it with the boundary conditions: 

3 2 

P = ^ ® ■ (2.8) 

yJ{drr]Y - 1 

where p = rji denotes a nucleation time. 

The tunnelling rate per unit four-volume is obtained 
from the imaginary part of the action with the complex 
ry, which becomes 


T{r]i) ~ exp[-2ImS'(7yj)]. 


(2.9) 


As a result [13|, the tunnelling probability in a spatially 
flat de Sitter universe is obtained as 


ImS^ = - 7 smm 




-ln(l + (3i7a/e)2) 


( 2 . 10 ) 


Note that the tunnelling rate is independent of the choice 
of T\i in the de Sitter spacetime. This expression reduces 
to the result for Minkowski spacetime known as the Cole¬ 
man’s bubble solution 0 ] 


ImS' = 


4e3 


in the limit of 77 —)> 0. 


B. Wall solution 

We discuss how the difference of the VEV’s (j)±, that 
is the difference of gravitational constants in the Jordan 
frame, affects the observational consequence. So we an¬ 
alyze the dependence of the ratio C?_|_/G_ = <1)^/$^ = 
exp[—-^(0_i_ — (j)-)] on the final result. For simplicity we 
discuss only the potential in the limit of e —0. That is 
the case of the potential for which both vacua have the 
same cosmological constant A 
The potential (12.51) in this limit is 

^o(</>) = ^(<(>"-m")' + A. (2.11) 


This gives 


G_ 


= exp 



( 2 . 12 ) 


So if we change the value of the parameter yx, we can 
discuss the dependence of the different gravitational con¬ 
stants. However, it also changes the height of the poten¬ 
tial barrier. In order to keep the same potential height, 
we adopt the following toy potential 

2A 

K)W = (-2 5)4 ('^~ + ]iy +A. (2.13) 


^ We give an analytic solution of a domain wall for small e in 
Appendix. 
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The values of the scalar field in the mother universe and 
our universe are given by 4>+ = bfi and (/)_ = —/i, respec¬ 
tively. The ratio of the gravitational constants is 


G_ 


= exp 


2(6 -f l)/r 


V7 


(2.14) 


Comparing (I2.11|) with the case 6 ^ 1 in (12.131) allows 
us to discuss the change of the ratio of G+ /G_ with the 
same height of the potential barrier (see Fig.l). Calculat¬ 
ing the wall solution and surface density of the wall cr for 
this toy potential and varying the value of 6(> —1), we 
can see the effect of the different gravitational constant. 

In order to evaluate the surface energy density a of 
a bubble, we consider a spherical symmetric and static 
solution (j){r). The basic equation is 


^ = 0. (2.15) 

r d<p 

Imposing the thin-wall approximation in which we can 
neglect the term of dr4>, we can integrate it as 


r dc^ 

4. \/2(Co(0) - Co(<(>z)) 


(2.16) 


where (pi = (t>(ri). 

For the potential (I2.13p . we can find the wall solution 
as 




with 


—/i -I- 




1 -I- tanh 


d = 


(1 + fc) 

2fi\/X 



(2.17) 

(2.18) 


Of course, in the limit of 6 = 1, we find that the solu¬ 
tion (12.171) is reduced to the usual wall solution with the 
potential (j2.11l) 


We can argue how the difference of gravitational con¬ 
stants in the Jordan frame affects on the bubble nucle- 
ation rate F, the wall dynamics and the primordial per¬ 
turbation through this surface density <7(6), which is a 
function of 6. 

Note that we find the same surface density a even when 
we take into account of small e-modification of the po¬ 
tential (see Appendix). 


III. DYNAMICS OF BUBBLE WALL 

In this section, we study dynamics of a bubble wall 
after nucleation of a bubble universe. To obtain an ana¬ 
lytic solution of a wall trajectory, we assume that a wall 
is infinitely thin and impose a junction condition across 
the wall. The bubble dynamics has been studied by sev¬ 
eral authors Sim, [l4 [ 2 ^. Following their works, we 
briefly summarize the results. 

Each spacetime is distinguished by the suffixes -I- and 
—, respectively. The bubble wall is represented as a time¬ 
like spherically symmetric hypersurface E dividing a false 
(-I-) vacuum and true (—) vacuum. We adopt a spa¬ 
tially flat slicing of de Sitter universe for both spacetimes, 
whose metrics are given by 

ds^ = —dt^ + exp(2i7±t±)(dr^ -|- r^dU^), (3.1) 

where H± = y^A±/ 3. Note that 

Hl-Hl= e/3. (3.2) 

When e 0, cosmological constants inside and outside 
of the wall have different values. Up to the first order of 
e, they can be approximated as A_ = A — e and A_|_ = A, 
respectively. The metric of the bubble wall takes the 
form 


, , r fiy/X, , 

(j)[r) = ^tanh —;^{r - n) 


(2.19) 


Then we obtain cr as 

pCO 

a = / dr 


/o 


-(9.0(r))2 + U(0) 


0+ 


d0v'2(Vb(/>)-Uo(0_)), (2.20) 


where we have used the thin-wall approximation. Sub¬ 
stituting (12.131) into the above equation, we obtain the 
result 


cr(6) = 4^(1+ 5). 


( 2 . 21 ) 


When 6 = 1, it gives the surface density of the original 
potential (12.111) as 


<7 6=1 = 


2^1 


We assume the stress-energy of this wall is given simply 

Sf2]y = , (3-4) 

where is the projection tensor on the 

hypersurface E, which is the metric tensor of (lT3ll . The 
unit spacelike vector normal to the hypersurface E is 
given by 

f dr dt \ f dr I dt 


From n^riu = 1, we obtain the relation equation 


_ 2 f _ 


dr 


\ dr J 


= 1 . 


(3.5) 


(3.6) 


By using x it can lead to 

dt± 1 


i-4(^) 


3 


( 2 . 22 ) 


(3.7) 
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In order to find the dynamical equation for the wall, 
we use the Israel’s Junction conditions as 


M;|e=0, (3.8) 

(3.9) 

where [X]!*! = X\+ — X\- for the variable X and S = 
The extrinsic curvature K^i, of the hypersurface E 
is defined by 


K^, = hP^hZV,np, (3.10) 

where Vo- denotes a covariant derivative with respect to 
g^y. The continuity equation is given as 

= (3-11) 

where || denotes a covariant derivative with respect to 
the projection tensor h^y. This equation is reduced to 
the dynamical equation for a as 


da{T) 

dr 


dr dt , . 


(3.12) 


where p and P denote the energy density and pressure 
in the bulk universes, respectively. If p + P 7 ^ 0, e.g. for 
the radiation or matter dominant universe, a changes in 
time. Throughout this paper, however, since we consider 
only de Sitter background universes, a turns out to be 
constant. 

From the junction condition (j3.8ll . we obtain the 
proper radius of the wall by 

R = a+r+ = a-T- . (3.13) 

To find the dynamical equation for R, we use the (90) 
component of the second junction equation ()3.9|) . which 
leads to 

K+-K^, = -^R\ (3.14) 

and its square becomes 

iKer = ^ [{K^ef - {Ktef - ■ ( 3 - 15 ) 

While K 00 is calculated as 

KI0 = |l + (3 jg) 

Using (I3.15|) . we can simplify the equation for R as 


where the second equality in (I3.18|) is found from Eq. 
(lO) . For c+, although e is assumed to be small, c+ can 
be positive because 

— 

V™PL 

can be smaller than unity, if we assume p <C mpL, where 
topl(= 1) is the reduced Planck mass. 

For de Sitter bulk universes, which we have assumed 
here, P is a constant and then a solution of a bubble 
radius R{t) is obtained by 

R{t) = — coshPr . (3.21) 

JD 

We can also solve the wall motion in the interior or 
exterior coordinates (t+, r_|_) or (t_, r_). In what follows, 
although we give the solution for the exterior coordinates, 
the same form of the solution is obtained for the interior 
coordinates just by replacement of + with —. 

From (I3.13p . we have 


4e 12 / 

^ “ (1 + 6)2 



Ini? = + lnr+, (3.22) 


and differentiate both sides of the equation with respect 
to T. Taking its square and using (13.6|) and (13.211) . we 
find the equation for (i(lnr+)/dr as 


d(lnr+) 9,(lni?)±i/2i?|c+| 

dr 1 - HlR^ ’ ^ 

which is integrated by using the solution of (13.211) as 


roo+ cosh Bt 
^ I sinh Bt ± c+ | ’ 


(3.24) 


where roo+ is an integration constant to be determined 
later and ± take +(—), if the sign of c+ takes a positive 
(negative) value. Inserting p.24l) into (I3.22p leads to the 
solution of as 


Hj^t+ = ln(| sinhPr ± c+|) — lnProo+ . (3.25) 

Finally we solve the wall motion as r+(t+), which de¬ 
notes a dynamics of the comoving bubble radius r for the 
observer in the exterior universe. 

Combining ^ ^ x -^{ar), (13.71) and (13.211) . we have 

a differential equation 

-P2j^2 ^ ^ ^ |c+|v/p2j^2 _ 1 _ (3 26) 

H 4 . r 


where 



B^R^{t) - 1 , 


(3.17) 


It is integrated when rewritten as equation for B^R? as 


r+(4) = Y/(a;i-l)2 + c;2i7-i. (3.27) 


c± 


P +(1 + 4 ) = P !(1 + c ^), 






(3.18) Here we take a nucleation time of a bubble as = 0, 
when we impose 9tr+(0) = 0. The solution takes the 

(3.19) 

same form for both interior coordinates too. 
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We can summarize dynamics of bubble wall in the in¬ 
side and outside observer as follows. At the initial time 
t± = 0, the bubble is created at the comoving radii 


The new potential U {ip) is assumed to be independent 
of V{(j)). In this case, we find the background equations 
as 


^’±( 0 ) = ri± 


e a 
^ ^ 4 


-1 


(3.28) 


and then expands. The comoving radii eventually con¬ 
verge to 



(4.2) 

n'-n^ = -\ip'\ 

(4.3) 

ip" + 2Hip’ + = 0 , 

(4.4) 


r±{t oo) = -I- c±ri± = roo± . (3.29) 

The integration constant roa+ given in (13.241) should be 
the same as the above value. The asymptotic comoving 
radii are rewritten by 

W = ^i?^" + (r.±)2. (3.30) 


where T-L = Ha and a prime denotes a derivative with 
respect to a conformal time 77 . We assume a slow-roll 
inflation for ip field for simplicity, i.e. 

a(t) = exp(i7t), H = --, -%-=H, (4.5) 

Tj 2 

where a dot denotes a derivative with respect to the cos¬ 
mic time t. We introduce slow-roll parameters as 


If we ignore small initial radii ri±, this equation basi¬ 
cally yields rao± ~ which means the final conver¬ 

gent radii are approximately described by the Hubble 
radii ^ in each coordinates. It is obviously noted that 
physical scale of the bubble radius continues to expand 
as A = a±r± oc , although the comoving radii con¬ 

verge to about the Hubble radii. 



£2 


g'l 

Hei ■ 


(4.6) 


Adding the perturbation of scalar field Sip and expand¬ 
ing the basic equations up to linear order, we find one 
master equation by use of the Mukhanov-Sasaki variable 
V, which is a linear combination oi Sip and tp. Defining 


IV. PERTURBATIONS 


M=^, V = a (^Sip + 


(4.7) 


In this section, in order to discuss observational effects 
of the different gravitational constant in the mother uni¬ 
verse, we analyze metric perturbations in the Newtonian 
(longitudinal) gauge. The perturbed metric is given by 


we find the basic equations for these variables as 

Au = z{v/zy, v = 9{u/ey, (4.8) 

where we have used 


= a^[—(I -I- 2'ip)drf' -|- (1 — 2'ip){dr'^ + r^dD^)]. (4.1) 

In the interior and exterior spacetimes, the tunnelling 
JBD scalar field takes the values of (p- and <()+, respec¬ 
tively, which provide non-zero cosmological constants. As 
a result, each bulk universes expand exponentially as de 
Sitter spacetime. However, such a de Sitter expansion 
will not end. In order to discuss more realistic cosmo¬ 
logical scenario, instead of a cosmological constant, we 
introduce another scalar field ip, which is confined at our 
vacuum state cp = (p-. This scalar field ip is responsible 
for inflation and will hnish the exponential expansion. 
The quantum fluctuation of this scalar field provides the 
density perturbations, which we will discuss here 


6 = 


n 

aip' 


1 _ I 

a^/2ei ’ 9 


(4.9) 


Then we obtain the closed master equation for u: 


u — 



(4.10) 


Under a slow-roll condition, the term 9” /9 is evaluated 
as 



2 1 £2 
" =4+^^ + y 


(4.11) 


For spherically symmetric perturbations, a solution 
u(r, 7?) = jr satisfies the equation 


^ In the outside mother universe, we may also find an infiaton, 
which will finish the exponential expansion. In our analysis of 
perturbations, however, we assume de Sitter expansion in both 
bulk universes. Hence our result does not change unless the 
outside inflation will end earlier than our inflation. There may 
be some effects on the perturbations from the infiaton field, which 
we ignore in our analysis. 


Am = (9^ -I- 2drlr)u = —k^u . 


The mode function with a comoving wave number k is 
obtained by 




U = 


\c^Hyy\-k7^) + c2Hyy\-kTi) 


(4.12) 
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where Ci ,2 are arbitrary constants. In the interior 
coordinate(—) system, v is reduced to Vk —>■ 
in the high-frequency limit, which corresponds to an adi¬ 
abatic vacuum. This condition fixes the above constants 
as (71,2. Then the solution u is determined up to a phase 
factor as 

(4,13) 

This mode solution is proportional to in 

the large |77_| limit, so interpreted as an incoming wave 
from the past. It propagates in the outward (large r_) 
direction from inside the bubble and then approaches to 
the boundary wall. 

This incoming wave is scattered at the wall and then 
divided into reflected and transmitted waves, which are 
written by 

Wrf = ^, (4.14) 

(4,15) 

respectively. Here the outside wavenumber: k+ is re¬ 
lated to as k+ = (a+/a_)fc_ because the phase factor 
^ik(r±r)) proper time dr^ = a^(d77^ — dr\) have 

to be invariant across the wall. 

The scattered mode Urf is moving in the inward direc¬ 
tion and then is going back to the interior bubble. While 
the mode Utr is transmitted to the outside bubble (the 
exterior universe). 

Expressing each perturbation as 

U- = Uin + /9urf, U+ = aWtr , (4-16) 

we match those wave solutions at the hypersurface by use 
of the junction condition. 

The matching has to be done for each proper time of 
the wall r. We then impose that the wave amplitudes 
and their normal derivatives are continuous at the wall 
hypersurface as 

M_(T) = u+(r), (4.17) 

ntd^u^{T) = n((ii9^u+(r). (4-18) 

The similar argument was given in as the matching 
perturbation across a hypersurface. 

From the wall motion (I3.24p and (I3.25P with (I3.27L 
the normal vector (1331) is found as 

j |c±| sinhHr — 1 ^ H^rooi coshiJr 

sinhHr -t- |c±| ’ 77±(sinhilr-f |c±|)^ 

(4.19) 


Rewriting a± and r]± as a function of r and using the 
matching conditions (14.171) and (14.181) . we obtain the am¬ 
plitudes of the reflected and transmitted waves. We plot 
some wave trajectory in Figl2] 

In the comoving radial coordinate r, the bubble radius 
(see the red line in Figl2]) increases from an initial value 
Ti to a convergent value roo± — 1/ H ±, while a comoving 
horizon 1 /{H±a±) (the green line) decreases monotoni- 
cally from an initial value 1 / H± where we have set = 1 
at an initial time t = 0 and approaches asymptotically 
to zero during de Sitter expansion. However if de Sitter 
expansion continues forever, a typical perturbation mode 
(as written in a blue line in FiglU) keeps to be located in 
superhorizon scale and cannot enter in the horizon again, 
which is not observable. Such a model does not describe 
a realistic scenario for cosmological perturbations. 


Hence as discussed above, we have added another 
scalar field </?, which is responsible for inflation. In a 
realistic inflationary scenario, inflation will end and lead 
to reheating of the Universe, connecting to a usual radi¬ 
ation dominated FLRW universe. Taking such a realistic 
scenario into account, we find the horizon scale increases 
again from the convergent value, and then such super¬ 
horizon perturbation will reenter into the horizon. The 
initial incoming wave of perturbations approaches to the 
wall, which is located on over the horizon (see the solid 
light-blue line with an arrow), and then is scattered into 
the reflected and transmitted waves. Now we shall ana¬ 
lyze how this reflection mode affects the observed quan¬ 
tities, i.e. the CMB power spectrum of the perturbations 
by the observer lived in the interior bubble universe. 


Generally speaking, if one of the bulk spacetimes is not 
de Sitter, it is difficult to find an analytic solution for the 
perturbations as well as the background spacetime. We 
need a complicated numerical analysis. However, since 
the universe expands exponentially during inflation, we 
assume the simplest situation for the calculation of the 
wave propagation, that is, the case of a de Sitter expan¬ 
sion in both sides of spacetime. 


In a de Sitter background, we can calculate the reflec¬ 
tion amplitude /3 and discuss how this mode affects on 
perturbations. We match perturbations in both sides at 
each proper time r. Eliminating a in the junction con¬ 
ditions (14.171) and (I4.18p . we find /3 as 

























0 Hr 1 



FIG. 2: We plot the trajectories representing the bubble wall (red) and horizon scale (green) as seen by a observer outside 
the wall in the units of H- = 1. The blue line denotes a typical perturbation scale, which is in superhorizon scale during de 
Sitter phase. It is shown by the light-blue lines that an incoming wave propagates from inside the bubble to the boundary 
wall, scattered at the wall and then divided into a reflected /3 and transmitted a waves. The pink line denotes incoming wave 
from the exterior bubble and scattered into a reflected /3 and transmitted d waves similar as light-blue lines. (Left) We plot 
them in the coordinate t and comoving radius r. (Right) We plot them in the coordinate conformal time r] = —l/{aH) and 
comoving radius r. It is shown that the slope of incoming wave (light-blue) is an angle of 45°. So at late time, the reflection is 
not generated since the slope of incoming asymptotically becomes a same slope of the wall trajectory, i.e., /3 ~ 0 and similarly 
all waves from the exterior bubble can be transmitted, i.e., fi ~ 1. 




-1 


+ 9+{T)H^Jli[-k+r]+)Hi^\-k_r]_) - g_{T)H‘;^l^{-k_r]_)Hl^\-k+ri+) 

X [(/+(t) - r_{T))H^J\-k+g+)H^J\-k_ig.) 

- 9+ {-k+r]+)Hl^'>{-k-r]- ) -b g- (T)iL,^+\(-fc_77_ {-k+g+) 

where * denotes a complex conjugate and 

/±(r) = (|c±| sinh^T - + (1 “ i^±’'oo±2/±(r)"^ coshSr)^ 2/±('r)"^ 

9±{t) = k±r^±{\c±\ sinh Bt - 1)?/^^(t), 
k±Broo± 


-fc±?7±(r) = 


H±y±iT) ’ 


with 


(4.20) 


(4.21) 

(4.22) 

(4.23) 


y±(r) = sinhSr + |c±| . 


(4.24) 


We note that in the above calculation we have used the same value of jz = 1/2 for the Hankel functions in both sides 
of the wall. Strictly speaking the order of the Hankel function has to be i'±, which depend on the slow-roll parameters 
ei, £ 2 . In our calculation, however, both sides of spacetime are approximated by de Sitter, that is, slow-roll parameters 
vanish and then we approximate the wave functions with v± = v = 1/2. We have also used the formula of the Hankel 
functions: dzH^{z) = vz~'^Hi,{z) — H^j^i{z). 

Using the definition of Hi/ 2 {z) and i? 3 / 2 (z), we find 


P = e 


2ik-r— 


f+{T) - f-jr) + -b i[g+{T) - g_(T)] 

_( f . (n-\ _ (n-\\ _ g+('^) I d-i'T 




+ - ^[9+{r) - g-(T)] 


{f+{r) - /l(r)) 


(4.25) 
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If e vanishes exactly, i.e., there is no energy difference between two vacuum states, since = iJ_, /+ = /_, 
5 + = g-, and 77 + = g-, then we find /3 = 0. However, we assume e/3 = ^ 0, which results in /3 oc e 0 if e 

is small. In fact, assuming e <C x (Tx/mpL)^ as well as e <C A, we hnd 


/3 


where we define 


g 2 *fe_ro ^ 


I \/l + ^0 
Co 


■ cosh(i?oT) ^I + - sinh(i?oT) M + 



(4.26) 


So 




.^tr 

co = ^, 


Co = 




Co 


vT+q 


% = 


(1 + cl)(. 


Co(sinh(Hor) + cq) ' 


(4.27) 


Similarly, we also have to estimate the contribution 
from the perturbation in the exterior universe (see the 
pink line in Figj^]). Now the incoming wave u is given as 

= (4.28) 

r+ 2 /c+ 

It is scattered at the wall and then divided into reflected 
and transmitted waves, which are written by 

= (4.29) 

r+ 2 A:+ 

utr = - -, (4.30) 


respectively. Expressing each perturbation as 

M+ = Itin + /3urf , U- = aUtr , (4-31) 

we match those wave solutions at the hypersurface by use 
of the junction condition. 

Eliminating /3 in the junction conditions (j4.17|) and 
(I4.18p . we find the transmitted waves of the perturbations 
in the exterior universe, that is moving inward direction 
to the interior universe. The transmitted rate a is given 

by 


77 + fc_r_ ,, 




V- k+r+ 


[f+ir) - /+(c)) {Hl^H-k+r]+)y 


X (/+(c) - fl{T)jHi^\-k+r]+)Hl:}'>{-k_r]_) - g+{T)Hl'^^^{-k+f]+)Hl:^'>{-k_r]_) 


1 -1 


Using the definition of Hi/ 2 {z) and 773 / 2 ( 2 ), we find 

^i7_ 


a = 


k+H+ 
When e is small, we find 


gi(fc_(r_+77_)-fc+(r++i7+)) ^ 


/+(r)-/;(r) 


(/+('c) - /-(c)) + + i[g+{T) - g-{T)] 


~ 1 


2ik-go H - (sinh(HoT) - Cq ^) 

v/r+72cosh(HoT)l ° ^ 


\/l + cp 


Co 


“ cosh(Hor) ( 4 - 3 + ^ ) - sinh(i?oT) (1 + 3 + ^') + ^ 


(4.32) 


(4.33) 


(4.34) 


Once we would obtain the solution of m, we have to convert it to the curvature perturbation and es- 






































Hence the curvature perturbation is given by 
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tiinate its power spectrum. In Fourier space, we find 

Cfc = "^klz = e^iukjoy. 


J 




dr 


k-v-^3/2i y2^77_ ^ 

dr 


Using the definition of Hankel function, a dimensionless power spectrum can be obtained by 


87r2 


TJ f n \ TJ 

{i + ik_7j_)—^p+ 1 + ik-Tj-^^ —— 


i/+ da 


77 _ i7_ d/3 77 + 

T dr ^£ 1 + dr 


(4.35) 


(4.36) 


where we neglected a slow-roll correction, i.e., £ 2 ± = 0 . 

We can divide the deviation from the perturbations in a standard inflationary model into three parts sq, si and S 2 , 
which are of order 0 {e), as 


Stt £ £i^_A^(fc_) = 1 -|- So -|- Si -l- S 2 4" Cl(£ ) 


(4.37) 


with 


So = 2Re 


—/3 + (5-l) 

£l_ 


Si = 2Re — - 


? 7 _ /ei_|_ df3 77 + da 


dr 


fllA y £i_ dr ^!Z± dr 


S2 = 27fc_77_ X Im I 


dr 

d'n_ 

d+ — 


V- 


dr)+ 

dr 


(5-1) 


(4.38) 


where Re and Im denote a real and imaginary part, respectively. 

At a superhorizon scale (fc_ 77 _ <C I), the contributions of so and si are important, which we can calculate 


So = 


ei'^coy/^ 


6 \/l + Cq cosh(HoT) 


sin( 2 fc ro) 


(sinh(RoT) - Cp 


+ 


\/TT^ 


Co 

— sinh(RoT) 


■ cosh(Ro''') 


cos(2fc ro) 






4- 


cos(2fc_ro) ^ 1 


and 


Si = 


e^^co7/o-v/eI+ j sin(2fc_ro) 


In— 2 + 
Co 


yo 


do 


2770 


770 cos(2fc_ro) 
coi 




(4.39) 


( 


cos( 2 fc ro) 


=) 


6\/l -I- Cp cosh(i?oT) I k-l^e\- co cosh^(i?oT) \/Ci- \/ei+ 

) ■ ^('“‘■(BoTjyiT^ - c„8mh(i3„r) + 1 + 


cosm(i?oT) 


1 + 


do' 


3r]o cosh(RpT)i/l + cj 


coyo^ 


\/£i+ 


(4.40) 
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where 


?/o = sinh(Bor) + cq . (4.41) 

The correction si is derived from the time derivative terms of j3 and a, which plays as the order 0(1/ cosh(i?o'r)) and 
then it can be ignored at the late time, although it can become a leading correction to (14.391) at the initial time. 

At the subhorizon scale (fc_77_>l), the dominant correction term is given by S 2 , which is evaluated by 


S2 = - 


e£^fc_?7oV^i+ j /cos(2A:_ro) 


6 


Co sin(2fc r'o) 




=v 


Co sinh(i?o'r) — 1 


V'T+cg cosh(Bo'r) 


Co 


\/Ci+^ fc-fya + cg cosh(i?oT) 

■ cosh(i?o'r)(^l + - sinh(i?o'r)(^l + 


2fc-7?o 
\/Cl + 

1 


?7o' 


Vo 

Co£ 


(4.42) 


In order to discuss how those correction terms can be 
observed, we shall depict them numerically with the spe¬ 
cific values of parameters, which we choose 

e=10-2£-^ andCT>0.4^-^ (4.43) 

This inequality guarantees to satisfy the condition 4e < 
3(7^ given as (13.201) . In FiglU we show a modulation 
factor of power spectrum (1 -f so + si + S 2 ) for different 
values of <t. It shows that the modulation factor can 
be larger in a subhorizon regime: k-^>l. For a large 
value of Co, the wall radius has still stayed to be constant. 
Namely cq > 1 gives a very tiny modulation, that is a 
nearly constant amplification in a power spectrum for 
T = I. Smaller value of cq shows that the modulation 
factor is increasing with oscillations in terms of k (e.g. 
see the plot for cq = 0.3), which amplitude gets larger 
for a smaller value of cq. 



k 




FIG. 4: All parameters except for co are same as used in Fig[3] 
(Top) We plot three contributions of so, si and S 2 for cq = 0.1. 
(Bottom) We compare the term 1 -I- S 2 with total one. The 
term S 2 can give us a good approximation for describing the 
total modulation function. The increasing oscillation appears 
around k-£ > 1. 


FIG. 3: We plot a modulation factor of power spectrum (1 -|- 
so -b Si -I- S 2 ) given in (14.371) . We set e = 0.01£~^ and t = £ 
for each value of co = ^cr/4 and assume slow-roll parameters 
on both sides of walls are equal as ^/si± =0.1. 

In FigHJ we plot how each correction contributes to 
the total modulation of power spectrum. We find that 
the term si can be ignored, while sq plays an important 


role in the super horizon scale k-£ < 1. On the other 
hand, S 2 plays an important role in the subhorizon scale, 
which gives an increasing function of k with oscillations 
for fc_£ > 1. This subhorizon effect comes from the one 
proportional to (fc_r 7 o)^ in (14.421) . That is transmitted 
wave propagating from the exterior universe. Note that 
the contribution of S 2 becomes larger than unity, which 
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means a breakdown of perturbative approach for e. The 
most effective term comes mainly from the one propor¬ 
tional to (fc_? 7 o)^ in (I4.42I1 . It gives the condition (I4.44|) . 
We can then evaluate a breakdown scale /cmax of our per¬ 
turbative treatment for e as 

< 0(1) -)■ /Cmax = 

(4.44) 

For Co = 0.1, we find a limited maximum value of 
wavenumber as fcmax — beyond which we cannot 

use this formula. 

Of course, these plots also depend on the parameter of 
the proper time r when the waves are scattered. However 
the plot shows a similar behavior when t < i, although 
we set T = i in Figl3] For cq = 0.1, the parameter t > 4£ 
gives no modulation of power spectrum since the wall 
radius converges to a constant value. 

In the observational point of view, the observed power 
spectrum would be generated by a wave scattered at the 
wall around which an e-folding number takes about 50 — 
60. If an inflationary period is much longer than 60 e- 
foldings, it maybe difficult to observe such effect on the 
CMB power spectrum because the observed modulation 
is generated at, t ^ i when the wall radius is almost 
constant. So we have to fine-tune the inflationary model 
whose period ends around 60 e-foldings. In this case, the 
value of the modulation is determined by the scattered 
waves at the wall position when inflation has just started 
and the radius of the wall is expanding. It will show the 
deviation from a standard inflationary model. 

Note that if the initial radius of the wall is much larger 
than the horizon scale H~^ = i, i.e., cq = cr£/4 <C 1, we 
may observe the original perturbation before the scatter¬ 
ing at large scale. So there may not appear the modula¬ 
tion for such a larger scale perturbation. 


ceyoij) 

V eP (1 -I- c§)^ ' 


V. THE EFFECT OF DIFFERENT 
GRAVITATIONAL CONSTANTS 


In the original Jordan frame, the gravitational con¬ 
stants are different in both vacua, which is given by the 
VEV’s of the scalar field $ by Eq. (12.21) . Since the scalar 
field (j) is given by Eq. (12.171) . the value changes from cj)- 
to (j)+ suddenly near the bubble wall. Hence the grav¬ 
itational constants are almost constant both inside and 
outside of the bubble except near the wall region. The 
VEV’s in a false(-l-) and true(—) vacua are written by 

(j)+=b^i, (/>_(= 0o) =-M- (5.1) 


Then we find 

$+ = 4>_ = 1, 


which gives 


G+ = 


(5.2) 


(5.3) 


Hence when we analyze the effect of the different grav¬ 
itational constant in a false vacuum, we should study the 
dependence of b on our results discussed in the previous 
section. If one wishes to set almost the same values for 
both gravitational constants, one needs to choose a neg¬ 
ative value of b such that 6 ~ —1, which implies that 
the value of </> in a false vacuum is close to that in the 
true vacuum ((^+). In Fig. [Sj we show the gravitational 
constants for two values of 6 {b = 1 and 2). 

1 


0.1 

log G 

0,01 


0,001 

■10 -5 0 5 10 

r-n 



FIG. 5: We plot the gravitational constants G for b — 1 and 
& = 2, by setting jj, — and A = 1. We have used the radial 
coordinate r in the Jordan frame, r; denotes the position of 
bubble wall. 


Since the parameter h is related to the surface density 
of the wall cr{b) as 


cr(6) 


3 


(1 + 5 ), 


(5.4) 


which is one of the key parameters in our results, we can 
discuss the effect of the different gravitational constant 
by changing a. In fact, the ratio of two gravitational 
constants is given by 


G- 


= exp 


12cr - 


(5.5) 


Now we study how change of the gravitational con¬ 
stant outside of our universe affects on the final results 
obtained in the Einstein frame in the previous section. 
The effects are found mainly in two points; one is the 
bubble nucleation rate and the other is the power spec¬ 
trum of curvature perturbations, which we shall discuss 
in due order. 


A. Nucleation rate of the bubble with 
different gravitational constant 


The bubble nucleation rate is given by (12.91) and (12.101) , 
which is described as 


F ~ exp 


TT^cr 


G. = Gn. 


2iJ3j ’ 


(5.6) 
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in the limit qof e —> 0. This result shows that the larger 
value of tr makes tunnelling process suppressed highly 
by a huge exponential factor. If cr ~ 0, i.e., 6 ~ — 1, 
which means the difference of gravitational constants is 
very small, the nucleation probability becomes very high. 
This is easily understood as follows: By assuming that 
/i and A are fixed, when the value of b changes from — 1 
to oo, the surface tension a takes the value from 0 to oo. 
It is because the wall barrier height is constant but the 
width becomes narrower as b decreases. When we change 
the gravitational constant in a false vacuum, however, if 
we fix cr, i.e., the potential is given by the relation such 
that 

- ct>_)\ (5.7) 

is constant, the nucleation probability does not change. 
Hence whether the universe with a different gravitational 
constant is plausible or not depends on the choice of the 
potential form. 


B. Curvature perturbations 

The power spectrum of curvature perturbations which 
reenter inside the horizon of our bubble universe is eval¬ 
uated by considering both reflected /3 and transmitted 
a waves. The modulation factor of a normalised power 
spectrum of curvature perturbations (|4.37I) - (I4.42I) is af¬ 
fected by changing cr(b). In what follows, we assume 
/i = just for simplicity, and then compare the case of 
6 = 1 with that of 6 = 2. If the surface density takes the 
value a = 0.4£“^ for 6 = 1, we find a = 0.6£~^ for 6 = 2. 
The gravitational constant G_|_ for 6 = 2 of the exterior 
of the bubble takes smaller value than the case of 6 = 1 
by the factor e~^ ~ 0.1 because (1 -I- 6) becomes three 
halves (see Figl5|). 

For small value of e, which we assume through this 
paper, which have been obtained as (I4.39|) - (j4.42[) . we 
plot the effect of difference of gravitational constants on 
the power spectrum in FiglHl The modulation factor 
(I-I-S 0 + S 1 +S 2 ), which describes a correction from a stan¬ 
dard scale-invariant spectrum, shows oscillatory bumps 
at a small scale k- >l/£ and its amplitude increases for 
smaller fc_. It also increases for a smaller value of the pa¬ 
rameter cr. By using (I4.42L the term such as cos(2fc_ro) 
gives the oscillation period Afc ~ n/ro oc a for small cr. 
Therefore when 6 becomes a half, the period Afc becomes 
more sharp by twice. 

We conclude the effect of different gravitational con¬ 
stant on the final power spectrum as follows. The typical 
signature appears on a small scale and it gives an in¬ 
creasing oscillation. Especially, if the gravitational con¬ 
stants for both sides of the wall is close to each other, i.e., 
cr Ri 0 (6 Ri — 1), the modulation of power spectrum be¬ 
comes larger. The feature where an increasing oscillation 
appears on small scale is generated by the transmitted 
d wave of perturbation in the exterior bubble universe 


(see (14.341) 1. and can be constrained by the observation 
of CMB. 



FIG. 6: We plot the final power spectrums for a = 0.4^ ^ and 
a = 0.6£~\ 


VI. CONCLUDING REMARKS 

We have argued a possible scenario that each universe 
takes its own different value of gravitational constant in 
the context of the emergence of the bubble universes. In 
order to construct a model that makes different gravita¬ 
tional constants corresponding to several vacua, we adopt 
the JBD theory, with which taking different values of the 
VEV of the JBD scalar field, the universes with different 
gravitational constants are achieved. In particular, we 
investigate the effect of different gravitational constant 
in the exterior universe on primordial curvature pertur¬ 
bations observed in the interior bubble universe. 

Via a conformal transformation, we can simply analyze 
the physics in the Einstein frame, that is the nucleation 
of bubbles, the evolution of bubble wall and matching 
curvature perturbations of the interior and exterior re¬ 
gion. In the Jordan frame, it becomes clear that the fi¬ 
nal power spectrum of curvature perturbation is affected 
by the bubble boundary through the surface energy den¬ 
sity cr of the wall. The key parameter is 6, which is 
6 = 0+/|^_|, where (/)± is the VEV’s of the JBD scalar 
field in the interior and exterior universes. The param¬ 
eter 6 is related to the ratio of gravitational constants 
as -2(6-1- l)/r/-\/7 = log(G+/G_). In Fig. O for differ¬ 
ent values of cr, we show the modulation of the power 
spectrum, which has several peaks with growing oscilla¬ 
tion in Fourier space. If the gravitational constant in the 
exterior universe gets smaller, 6 or cr increases and then 
the modulation becomes smaller. Choosing a smaller pa¬ 
rameter cr will heighten the amplitude of peak shown in 
the modulation of power spectrum of curvature pertur¬ 
bations. 

In our result, we assume that the e-foldings of inflation¬ 
ary period is about 60, because the initial modulation of 
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perturbations will not be observed for inflationary models 
with the e-foldings more than 60. In this paper, we have 
not discussed after inflation. The result obtained from 
perturbations matching at later time such as a radiation 
dominated era will be more interesting and important, 
although it is difficult to describe a trajectory of the bub¬ 
ble wall in an analytic form. In order to solve dynamics 
of wall in this phase, we may need numerical treatment. 
Along this direction, we hope to make a progress in our 
future work. 
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Appendix A: domain wall solution and surface 
energy for small e 


and expanding it up to the first order of small parameter 
e, we find the following equations: 


I ,, dVo 


f>0 + -(t>0 - = 0 ^ 

r d(p 


—^(0o)0i- 


d(j)^ 


(6-b I)/r 


By using thin wall approximation, Eq. 
0-th order domain wall solution as 


00 = —M ■ 


( 6 + 1 ) 


2 ^ 2 
Using this solution, we find 

^(0o) = T^-I- 


fi tanh 


2 d 


d(ir 


{h + iy 


2cosh^ (^) 


(AS) 
0. (A9) 
gives the 


(AIO) 


(All) 


If we ignore the deviation form unity near r = this 
term can be treated as a constant. Then we find an 
approximate solution for (j)i as 


( 6 + 1 ) . D 

4A/r3 [ dcosh^ (l^) 


(A12) 


where D is an arbitrary constant. 

This solution (j) = (pQ -\- e0i can be rewritten as 


We consider the following potential 


V = Vo + eVi 
2A 


( 6 + 1 ) 

e 

+ - 


(6 + l)/i 
There are three extrema 


4 (0 - bfj.) (0 -b fj.y 
(0 — b^) + A . 


(Al) 


0-k 

, (6+1) 

= ‘ 4Ap3 ' 

(A2) 

4>b 

6-1 (6+1) 

- 2 ^ + " 2AA.3 ’ 

(A3) 

0- 

(6+1) 

= ^ ^ 4Am3 ■ 

(A4) 


0+ and 0_ correspond to the local minima, while 0 b is 
the local maximum, at which the potential values are by 
given 


Although (j)+ can be shifted by changing a free parameter 
b, the potential barrier and the minimum values do not 
depend on b. 

Now we hnd a domain wall solution for this potential 
V. Plugging (p = 00 (r) + e0i(r) into the basic equation 


(5-1) (6+1) 

-/i e- 


2 

( 6 + 1 ) 


4A//3 


/rtanh 


dr. 


2 d 


(A13) 


where = -j^D. This domain wall solution is just 
shifted from one with e = 0 by and —Sr^ in p- 

and r-directions, respectively. The domain wall structure 
itself does not change. 

Next we shall calculate the surface energy density, 
which is defined by 


dr 


-(drpy + Udw(0) 


(A14) 


where Vdw is the part of the potential which contributes 
to the structure of a domain wall. 

Plugging the above solution (IAI3|) into the kinetic term 
and potential, we find up to the first order of e as 


I 


{drpy = 


Xfjp 



= A + 0(e2), 

(A5) 

z 

v{4>b) 

= A + - 2 A ^(^^) ’ 

(A6) 

U(0) 

U(0_) 

= A-e + 0(e2). 

(A7) 



8 cosh 


4 / r—rj+<5re 


(A15) 


2d 


X^jT 


+anh(%f-) 


+ 


8 cosh‘ 
e 


4 r—ri~\-5r, 

tanh 


2d 

r - Ti 


4cosh^ (^) 


2 d 


- I 


+ A. 


(A16) 


The potential V contains a back ground bulk energy, i.e., 
cosmological constants; 


Wuik = A - - [I - sgn(r - Ti)] . 


(A17) 
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where the sign function sgn(a;) is defined by 

/ / 1 if a; > 0 

sgn(x)= if ^<0 . 


(A18) 


Since the bulk energy should not be included in the sur¬ 
face energy of the wall, we define 


kow := — khuik 


_ etanh(^) 

8 cosh^ (4 cosh^ {^) 


tanh 


2d 
r -n 
2 d 


- sgn(r - n) 


(A19) 


In the thin-wall approximation (d <C Tj), the surface en¬ 


ergy (T is evaluated approximately as 


cr 



-{dr(j)Y + ^Dw(<^) 


dr 

~ ~T J_^ cosh" ’ 

which gives the same result as the case with e = 0, 


(A20) 

i.e., 


^/Xp?{b + 1 ) 


(A21) 


It is plausible because the domain wall structure does not 
change by the small e-modification of the potential. 
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